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Abstract The paper examines a trapped one-dimensional system of multicomponent spinless fermions 
that interact with a zero-range two-body potential. We show that when the repulsion between particles 
is very large the system can be approached analytically. To illustrate this analytical approach we 
consider a simple system of three distinguishable particles, which can be addressed experimentally. For 
this system we show that for infinite repulsion the energy spectrum is six fold degenerate. We also 
show that this degeneracy is partially lifted for finitely large repulsion for which we find and describe 
corresponding wave functions. 

Keywords strongly interacting few-body systems • one dimensional harmonic traps • multicomponent 
fermions ■ Tonks-Girardeau gas 


1 Introduction 

Experimental study of few-body physics in cold atomic gases is a complicated task since such systems 
usually have a relatively large particle density and, hence, many-body correlations should be taken into 
account. Only very recently setups with small particle numbers were realized in Heidelberg Qii 
where ground state systems of a few fermionic atoms, were prepared in a quasi-one-dimensional 
trap. Such setups pave the way for the experimental study of few-body correlations where accurate 
theoretical description can be obtained through advanced numerical investigation ia. This paper 
overviews the newly developed method Q to study the mentioned experimental setups in the limit of 
strong repulsion between particles without applying coimlicated numerical routines. This analytical 
approach gives a description of the Tonks-Girardeau gas aSQ of a few particles in a trap. Moreover, 
as a by-product results obtained with this method can be used as a reference point for numerical 
calculations. 

The structure of the paper is the following: in section l^we introduce the Hamiltonian that is widely 
used to describe the relevant experimental setups BiiiiKini , in section [3] we illustrate the 
approach using the simple system of three distinguishable particles in a harmonic trap, that to the 
best of our knowledge was not addressed before. There we find and describe eigenstates of such system 
in the limit of strong interparticle interaction. 
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2 Formulation of the problem 


This paper considers N particles of equal mass m in one spatial dimension. Additionally, the following 
assumptions are applied: i) the particles can be divided into classes of identical spinless fermions; ii) 
the system is trapped by some external potential, Vext] Hi) the interparticle interaction is assumed to 
be of zero range, V = gS{xi — Xj), where g is a strength parameter and Xi^xj are the coordinates of 
particles i and j. These assumptions lead to the following Hamiltonian 
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This model Hamiltonian is often used to stut^ the relevant experimental setup of a few fermionic atoms 
in different hyperfine states This model Hamiltonian allows one to explore 

analytical approach only in limiting cases, e.g. two particles in a harmonic trap [l^ or N particles 
without an external confinement |13| . such that the theoretical study for the relevant experiments is 
usually provided using different numerical techniques, e.g. in. It is shown in Q that one can find 
eigenstates for such Hamiltonian in the case of very strong repulsion between particles, i.e. if I /5 —>■ 0. 
This regime is experimentally accessible Q and very interesting from the theoretical point of view 
since it allows a theoretician to gain knowledge about strongly interacting systems. 

In this paper we overview the method from Q using a simple system of three distinguishable 
particles trapped in a harmonic oscillator trap. In other words we consider the Schrodinger equation 

(E ^ E + d^ixi - X2) + gS(x2 - X3) + gS(xi - xsjj <1' = EW, ( 2 ) 


where E is the wave function, E is the energy of the system, uj is the frequency of the harmonic 
oscillator, g is assumed to be large and positive, such that 1/g ^ 0. Eq. ([2]) can be recast into the free 
Schrodinger equation plus the boundary conditions 
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where i,j = 1,2,3 and i ^ j. Eqs. (I2|) and (I3|) contain all ingredients that are needed to describe a 
general method moreover, the system of three distinguishable particles is interesting on its own 
rights as experimentally relevant. 


3 Method 

3.1 Case of infinite repulsion between particles, 1/g = 0. 

To solve eq. (I2|) for large values of g we first notice that if 1/g = 0, then eq. ([31) yields = 0, 

which means that eq. ([3|) can be seen as an equation for free particles that cannot penetrate through 
one another. This problem formally resembles the case of three identical spinless fermions with the only 
difference that now the total wave function not need to have a continuous derivative when two particles 
meet. This possibility to have a discontinuous derivative follows from eq. However, since the 
particles cannot penetrate through one another it is enough to solve the Schrodinger equation only for 
a given ordering of the particles, i.e. xi < X 2 < X 3 , with the boundary conditions 'i'xi=x 2 =X 3 — 0. 

These boundary conditions can be satisfied only for the energies from the eigenspectrum of spinless 
fermions {Eq, Ei,...}, which implies that the total wave function for 1/^ = 0 can be built using the 
wave function of spinless fermions Ep for a given energy Ei 

' aiEp for X 2 < xi < X 3 
a 2 'I'F for X 2 < X 3 < xi 
a 3 <I'F for X 3 < X 2 < xi 
a^Ep for X 3 < xi < X 2 
a^Ep for xi < X 3 < X 2 
^ uqEf for xi < X 2 < X 3 


E = { 


( 4 ) 
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where are real coefficients. To obtain the wave function 0 we assumed that the energy Ei corre¬ 
sponds to only one wave function Ep. The extension of the method for the more general situation when 
the energy Ei corresponds to more than one wave function is discussed in refs. mM- The possibility 
to build the wave function in the form of eq. (|3]) implies that for 1/g = 0 the energy spectrum of eq. ([2]) 
is six fold degenerate, since all six coefficients Oi in eq. 0 are linearly independent. 


3.2 Case of finitely large repulsion between particles, 1/g ^ 0. 


Method. To find the behavior of the energy in the vicinity of 1/g = 0 we use the Hellmann-Feynman 
theorem to obtain the derivative of the energy 

dE ^ J dxidx2{'I')l^=^^ + f dxidx3{E)l^^^^ + f dx2dx3{E)l^^^^ 

dg {E\E) ■ 

Next we notice that from eq. ([3]) it follows that 'i'\xi=xj 1/g + o(l/g), which allows us to conclude 
that E Ei — K/g o(l/g), where the parameter K is given by 
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The limit in this equation can be easily taken, since {'I')xi=xj ~ 1/5^ + oil/g"^), which yields 
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Using the wave function (jj]) the value of K can be written as 


K = -! 
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where 7 is defined as 
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To determine the behavior of the energy we need to find the coefficient K. To do so we use the 
variational treatment, i.e. we minimize the energy by varying K with respect to the coefficients a^. 



( 10 ) 


In this way we find the wave functions (defined by Oi and eq. (|3])) to which the states outside of the 
degenerate point (1/g = 0) are adiabatically connected. This procedure is variational, since it relies on 
varying K with respect to the coefficients a^. However, it produces an exact solution, since the wave 
function E for 1 /g = 0 can always be written in the form (jT]), so we vary in the full space. 

Solution. Equation m produces the system of linear equations 


-06 -I- 2ai — 02 = Kai/'y , 

—Oi-i -I- 2ai — Oi+i = Koi/j, for j = 2, 3,4, 5 
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This can be seen as a usual eigenvalue problem for a system of linear equations, which can be easily 
solved using standard techniques from linear algebra. The solution contains six eigenstates: 

= 47 —>■ 02 = —Oi, 03 = Oi, 04 = —Oi, 05 = Oi, 06 = —Oi 

^^2 = 37 —>■ 02 = —ai/2,03 = —ai/2, 04 = oi, 05 = —oi/2, oe = —oi/2 
= 37 —>■ 02 = —2ai, 03 = 01,04 = 01,05 = —2ai, oe = oi 
= 7 —^ ^2 — ai, 03 = 0,04 = —Oi, 05 = —Oi, 06 = 0 
=7 —> 02 = 0,03 = —oi, 04 = —oi, 05 = 0, 06 = oi 
Kq = 0 —y 02 = 01,03 = 01,04 = 01,05 = oi, 06 = oi 

where the coefficient oi is left for the overall normalization. First of all one notices that the lowest 
energy state defined by Ki has a fully symmetric wave function and can be obtained from the Fermi- 
Bose mapping, which is first discussed in ref. This is always the case for the Hamiltonians in the 
form of eq. o, without any symmetry requirements. The state with Kq = 0 corresponds to a fully 
antisymmetric wave function of the system of spinless fermions. Four other states are the same as for 
the system of two spinless fermions and a third particle Bii- The degeneracy of the spectrum, 
K 2 = K 2 , and K 4 = K^, comes from the observation that one can antisymmetrize either xi and X 2 
particles or xi and 0:3 particles. Combinations of states with K 2 , K 3 , K 4 , also appear for two spinless 
bosons and a third particle as shown in [T^. 

In this way we obtained the total energy spectrum of the system of three strongly interacting par¬ 
ticles, where six states can be divided into three classes: three spinless bosons, three spinless fermions, 
and two spinless fermions interacting with a third particle. 


4 Conclusions 

In this paper one-dimensional strongly interacting systems of particles were investigated. We presented 
the method to treat such systems with zero-range repulsive interactions of large strength. To illustrate 
the approach we analyzed a simple system of three particles. For such a system we found the energy 
spectrum and the corresponding wave functions. 
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